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ABSTRACT 

Standing slow mode waves in hot flaring loops are exclusively observed in spectrometers and 
are used to diagnose the magnetic field strength and temperature of the loop structure. Due 
to the lack of spatial information, the longitudinal mode cannot be effectively identified. In 
this study, we simulate standing slow mode waves in flaring loops and compare the synthesized 
line emission properties with SUMER spectrographic and SDO/AIA imaging observations. We 
find that the emission intensity and line width oscillations are a quarter period out of phase 
with Doppler shift velocity both in time and spatial domain, which can be used to identify a 
standing slow mode wave from spectroscopic observations. However, the longitudinal overtones 
could be only measured with the assistance of imagers. We find emission intensity asymmetry 
in the positive and negative modulations, this is because the contribution function pertaining to 
the atomic emission process responds differently to positive and negative temperature variations. 
One may detect half periodicity close to the loop apex, where emission intensity modulation is 
relatively small. The line-of-sight projection affects the observation of Doppler shift significantly. 
A more accurate estimate of the amplitude of velocity perturbation is obtained by de-projecting 
the Doppler shift by a factor of 1 — 26>/7r rather than the traditionally used cosO. If a loop 
is heated to the hotter wing, the intensity modulation could be overwhelmed by 
background emission, while the Doppler shift velocity could still be detected to a 
certain extent. 
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1. Introduction 

Magnetohydrodynamic (MHD) waves are be¬ 
lieved to play a significant role in the formation 
and dynamics of the solar atmosphere. They may 
contribute significantly to coronal heating (see re¬ 
views by Klimchuk 2006; Taroyan & Erdelyi 2009; 
Parnell & De Moortel 2012; Arregui 2015) and 
solar wind acceleration (see e.g., Ofman 2010; 
van der Holst et al. 2014). During the past 
decade, a number of MHD wave modes of coro¬ 
nal loops were detected with modern instruments, 
e.g., standing and propagating fast kink mode 
(Nakariakov et al. 1999; Aschwanden et al. 1999; 
Williams et al. 2002; Van Doorsselaere et al. 2008), 
fast sausage modes (Asai et al. 2001; Melnikov 
et al. 2005), standing and propagating slow mode 
(Wang et al. 2003b,a; De Moortel et al. 2000, 
2002b,a; Wang et al. 2009a,b; Yuan & Nakariakov 
2012; Krishna Prasad et al. 2014). 

MHD wave theory in structured plasma forms 
a solid basis for a wave-based plasma diagnostic 
technique - MHD coronal seismology (see Nakari¬ 
akov & Verwichte 2005; De Moortel & Nakariakov 
2012, for recent reviews). MHD seismology was 
successfully applied in estimating the coronal mag¬ 
netic field (Nakariakov & Ofman 2001), transverse 
loop structuring (Aschwanden et al. 2003), Alfven 
transit times (Arregui et al. 2007), polytropic in¬ 
dex (Van Doorsselaere et al. 2011b), thermal con¬ 
duction coefficient (Van Doorsselaere et al. 2011b), 
the magnetic topology of sunspots (Yuan et al. 
2014b,a; Jess et al. 2013), and the magnetic struc¬ 
ture of large-scale streamers (Chen et al. 2010, 
2011). It could also be used to deteruiiue the 
coroual deusity scale height (Audries et al. 

2005), to quautify the expausiou factor of 
the coroual loops (Verth et al. 2008), aud 
to probe the characteristic spatial scale of 
raudouily structured plasuias (Yuau et al. 
2015). 

The slow sausage mode was initially theorised 
by Edwin & Roberts (1983). It is a compressive 
mode characterised by axisymmetric longitudinal 
displacement of the plasma fluid. Gas pressure is 
the main restoring force. In a \ow-/3 plasma, this 
mode does not cause significant contraction or ex¬ 
pansion of the loop cross-section, nor a displace¬ 
ment of the loop axis. Standing slow modes are 
frequently observed as oscillations in the plasma 


emission intensity and Doppler shift velocity in 
hot flaring coronal loops (> 6MK, see review by 
Wang 2011). They are exclusively detected in 
the hot emission lines, i.e., Ee XIX and Ee XXI 
lines as recorded by the Solar Ultraviolet Measure¬ 
ments of Emitted Radiation (SUMER) spectro¬ 
graph onboard Solar and Heliospheric Observatory 
(SOHO) (Wang et al. 2002, 2003b, a), and S XV 
and Ca XIX lines observed by the Bragg Crys¬ 
tal Spectrometer (BCS) onboard Yohkoh (Mariska 

2006) . Only recently, the Solar Dynamics Ob¬ 
servatory (SDO)/Atmospheric Imaging Assembly 
(AIA) imager observed, in the 13lA bandpass, 
multiple reflections of a propagating slow wave in a 
hot coronal loop (Kumar et al. 2013). A slow com¬ 
pressive mode is found to be launched by repeti¬ 
tive magnetic reconnections occurring at one of 
the footpoints (Kumar et al. 2015). 

Quasi-periodic pulsations (QPP) in solar and 
stellar flares are thought to be caused by MHD 
waves (Nakariakov & Melnikov 2009; Anfinogen- 
tov et al. 2013). Oscillations with periods at tens 
of minutes are ascribed to modulations by slow 
mode MHD waves (Van Doorsselaere et al. 2011a), 
while short period (sub-minute) oscillations are 
suggested to be modulated by fast mode waves 
(Van Doorsselaere et al. 2011a; Kupriyanova et al. 
2013). Recently, Kim et al. (2012) detected oscil¬ 
lations during a Ml.6 flare with the Nobeyama Ra- 
dioheliograph (NoRH) and the 335 A EUV chan¬ 
nel of SDO/AIA. The intensities in radio and EUV 
emissions oscillates with a period of about 13.8 
minutes and are damped with a decay time of 
about 25 minutes. Kim et al. (2012) inter¬ 
preted the detected ten-minute time scale 
oscillations as standing slow magnetoacous¬ 
tic waves. SUMER-like oscillations, but 
with shorter periodicities, are detected in 
soft X-ray emissions (Ning 2014). 

The standing slow mode wave oscillates with 
a period of about ten minutes and with a veloc¬ 
ity amplitude of a few tens of kilometers per sec¬ 
ond (Wang et al. 2003a; Wang 2011). Significant 
oscillation is normally detected at the loop apex 
and becomes absent at the footpoints (Wang et al. 

2007) . This is consistent with a scenario that an 
anti-node in the density perturbation of the funda¬ 
mental standing slow mode is located at the loop 
apex (Wang et al. 2007). This kind of waves are 
damped within a few oscillation cycles, which are 
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believed to be caused mainly by thermal conduc¬ 
tion in hot flaring loops (Ofman & Wang 2002; 
Selwa et ah 2005). 

The fundamental slow standing mode 
of a hot loop appears to be triggered by 
asymmetric heating at one footpoint and 
is rapidly established within one wave cy¬ 
cle (Wang et al. 2005). Pressure pulses 
launched close to a footpoint cannot excite a 
fundamental slow mode fast enough to com¬ 
pensate the strong damping (Selwa et al. 
2005). Taroyan et al. (2005) showed ana¬ 
lytically that the fundamental slow mode 
of a 6MK-loop could rapidly be excited 
by a single impulsive heating with a time 
scale that matches the loop period. For¬ 
ward modelling using simple one dimen¬ 
sional hydrodynamic model were performed 
to distinguish the propagating and stand¬ 
ing slow waves in both cool and hot loops 
(Taroyan et al. 2007; Taroyan Bradshaw 
2008). They also reported that Doppler 
shift variation is a more reliable observable 
to detect a slow wave, while intensity mod¬ 
ulation would be phase-shifted by heating 
(Taroyan et al. 2007) or contaminated by 
the background plasma emissions (Taroyan 
Sz Bradshaw 2008). 

MHD seismology with the standing slow mode 
was applied successfully to estimate the magnetic 
field strength and the time-dependent plasma tem¬ 
perature of a coronal loop (Wang et al. 2007). 
The application of MHD seismology relies on sev¬ 
eral nontrivial factors, e.g., the analytical model, 
mode identification, line-of-sight (LOS) effect, 
and plasma emission. Forward modelling was 
attempted to help interpreting the observations 
correctly. Cooper et al. (2003a,b) investigated the 
LOS effect on imaging observation of the emission 
intensity variation of fast kink and sausage mode, 
and explained the intensity perturbations of fast 
wave trains observed by the Solar Eclipse Coronal 
Imaging System (SECIS) instrument (Williams 
et al. 2001, 2002). Gruszecki et al. (2012) per¬ 
formed a three dimensional numerical simulation 
of the fast sausage mode of a plasma cylinder 
and investigated the geometric effect of simple 
LOS integration. Antolin & Van Doorsselaere 
(2013) and Antolin et al. (2014) developed ad¬ 
vanced models and included atomic emission ef¬ 


fects, using the CHIANTI atomic database (Dere 
et al. 1997). They found that the LOS effect 
and spatial resolution could significantly affect 
the intensity modulation and spectral character¬ 
istics of the fast sausage mode. Reznikova et al. 
(2014, 2015) used the same model and further in¬ 
vestigated the gyrosynchrotron emission intensity 
variation using the East Gyrosynchrotron Codes 
(Eleishman & Kuznetsov 2010). The radio emis¬ 
sion intensity of the fast sausage mode oscillate 
in phase for all frequencies, while for certain LOS 
angles, the optically thick and thin radio emis¬ 
sion are anti-phase. It opposes previous find¬ 
ings of Mossessian & Eleishman (2012), which did 
not consider the inhomogeneity of the emitting 
source along the LOS. Kuznetsov et al. (2015) 
used a semi-torus model to forward model 
the gyrosynchrotron radio emission of both 
propagating and standing slow modes in a 
curved magnetic structure. 

In this study, we perform forward modelling of 
standing slow modes of hot flaring coronal loops 
and predict their spectroscopic and imaging ob¬ 
servational signatures. We use the slow wave 
model in hot coronal loops (> 6 MK, see Section 2) 
and utilize the CHIANTI v7.1 atomic database 
(Dere et al. 1997; Landi et al. 2013) to synthe¬ 
size plasma emission in the SUMER Ee XIX line 
and SDO/AIA 094 A bandpass (see Section 3). 
Then the results and conclusions are summarised 
in Section 4 and Section 5, respectively. 


2. Model 


2.1. Standing slow mode 

In this study, we only consider a standing slow 
mode in a simple plasma cylinder embedded in a 
uniform plasma. The magnetic field is par¬ 
allel to the axis of the plasma cylinder (i.e., 
z-axis), Bo = Bqz. The equilibrium mag¬ 
netic field Bq, plasma density po and tem¬ 
perature To are the piecewise functions of 
the r-axis: 


Bo, po^ To 


Bi,pi,Ti :r<a 

Bq , Pe: T0 . T ^ Qj, 


where a is the radius of the loop. Hereafter, 
we use subscript T’ and ‘e’ to differentiate 
the internal and external equilibrium values 
of the loop system. Effects of plasma strati- 


3 


fication and loop curvature are ignored. We focus 
on observational features caused by optically thin 
plasma emission, LOS integration and instrument 
response function. We limit our study to the ax- 
isymmetric mode m = 0 (sausage mode, see Edwin 
& Roberts 1983). 

The linearised ideal MHD equations (see, e.g., 
Ruderman & Erdelyi 2009) give the perturbed 
quantities on top of the magnetostatic equilibrium: 

Pi = -v (2) 

= + — [(Bo • V)bi + (bi • V)Bo], 

Ot jlQ 

( 3 ) 

bi = V X (C X Bo), (4) 

Pi - Clpi = ^ ■ (QWpo - Vpo), (5) 

where ^ is the Lagrangian displacement vector, 

Po is the equilibrium plasma pressure, pi, pi and 
bi are the perturbed plasma density, pressure 
and magnetic field, Pti = Pi + bi • Bq/po is 
the perturbed total pressure, po is the magnetic 
permeability in free space. We define the key 
characteristic speeds to describe the loop sys- 
tem, Cs = i/lPo/po, Ca = Ba/^papa, Ct = 
CaC.IsJCI + Cl are the acoustic, Alfven, and 
tube speed, respectively (Edwin & Roberts 1983); 
cjs = uja = CAk, ujt = CTk are the corre¬ 
sponding acoustic, Alfven, and tube frequencies, 
respectively, where k = nnlLo is the longitudinal 
wavenumber, n is the longitudinal mode number, 

I/O is the length of the loop, 7 = 5/3 is the adia¬ 
batic index. 

Equation 2-5 are solved in cylindrical coor¬ 
dinates (r, 0 , z) with the boundary condition at 
r = a where the radial displacement and 
the total pressure are kept in balance. In case 
of the standing slow mode with m = 0 , we 
Eourier-analyse the perturbed quantities by as¬ 
suming Pti = A7l{r) cos{ujt) cos{kz), where A is 
the amplitude of the perturbed total pressure, and 
IZ is defined as 

7 ? - ^ 7 ?' (fW 

PoiCi + Cl){LO^-co^^)dt- 

Here 7Z' is the relevant ‘P’ in Antolin & Van 
Doorsselaere (2013). We use a cos{kz) longitu¬ 
dinal profile in total pressure perturbation, so 
the density (temperature) nodes are fixed at foot- 


points, while the longitudinal velocity perturba¬ 
tion follows a profile of sm{kz) and thus has a 
node at the loop apex for the fundamental mode. 
A key derivation is that 




cj^Pti 


poici + cDioj^ - co^y 


( 7 ) 


The perturbed total pressure must satisfy 


d^-Pri 

d ^2 


dPri 

rdr 


— k^Pti 


= 0 , 


( 8 ) 


where ^ 2 ^ ^^ is a modified ra- 

dial wavenumber and has the dimensionality of 
wavenumber k. Equation 8 holds for both inter¬ 
nal and external plasmas, where all quantities are 
piecewise functions of r (Figure 1). Equation 8 
gives 


f : r ^ a 

\ Ko{Krer) :r > a. 


(9) 


By matching the boundary conditions, we ob¬ 
tain the dispersion relation for the fast and slow 
sausage body mode (Edwin & Roberts (1983), 
> 0 and < 0, hence we re-define \Kri\ = 


Hvre I^o(^'7’e^) ^ri ^od^T’il^) 

pei^Ae ~ Ko{hireO) Pi{^Ai “ ’ 

( 10 ) 

where Jq is the 0-th order Bessel function of 
the first kind and Kq is the 0-th order modified 
Bessel functions of the second kind; Jq and are 
the corresponding derivatives with respect to /^^r. 
The perturbed thermodynamic quantities, which 
affect plasma emissions, are the velocity vector 
V = d^/dt, the plasma density p, and the temper¬ 
ature T: 


Vr — — 


AdlZjdruj 
Pq{lO^ -Lol) 

ATZC^kco 
poVliuj'^ - w|) 


sin(ci;t) cos{kz), 
sin(co’t) sm{kz)^ 


v<f) = 0 , 


( 11 ) 

( 12 ) 

(13) 


Pi = 
Ti = 


AIZ uj‘^ 

An{-f - l)To 

Io(C2 + Ci) (u;2 


cos(co’t) COs(/C2(), (14) 

CJrp) 

^2 

-^ cos(co’t) cos{kz), 

— CJrp) 


(15) 
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where Tq is the equilibrium plasma temperature. 
The thermodynamic variables are related by the 
equation of state for fully ionized hydrogen p = 
2kBpTfmp, where p = po + pi, T = Tq +Ti, /cb is 
the Boltzmann constant, rup is proton mass. 

Equation 8-15 holds for both the fast and slow 
sausage mode. The solution depends on whether 
the equations are solved in the Alfven [c^Ai^^Ae] 
or acoustic [cjxi^^si] frequency range (Saku- 
rai et al. 1991). We note that Vr/vz oc — 
— ujjJ. Therefore, in the Alfven fre¬ 
quency range, Vr/vz ^ 1 , it corresponds to the 
fast sausage mode (Antolin & Van Doorsselaere 
2013; Reznikova et al. 2014); while in the acoustic 
frequency range, Vr/vz <C 1 , this corresponds 
to the slow sausage mode (see, e.g., Moreels & 
Van Doorsselaere 2013), which is the topic of this 
study. 

2.2. Hot flaring coronal loop 

Hot coronal loops are complex and highly 
dynamic structures heated to 2-30 MK by 
flares (see a review by Reale 2014). A coro¬ 
nal loop may have unresolved fine struc¬ 
tures (Priest et al. 2002; Van Doorsse¬ 
laere et al. 2014) or multiple-strands (Pe¬ 
ter et al. 2013; Antolin Rouppe van 
der Voort 2012; Antolin et al. 2014, 2015; 
Scullion et al. 2014). Heating/cooling (see 
e.g., Klimchuk 2006; Hood et al. 2009) and 
the associated flows (see Winebarger et al. 
2002) are usually detected at the footpoints 
of the loops. MHD wave theory (see e.g., 
Edwin Roberts 1983; Sakurai et al. 1991; 
Goossens et al. 2011) normally assumes 
that a quiescent loop is in equilibrium with 
the ambient plasma, therefore the heating 
and cooling time scale should be sufficiently 
longer than the MHD time scales (in order 
of minutes in the case of slow modes). If 
a loop is subject to active heating or cool¬ 
ing, the Wentzel-Kramers-Brillouin (WKB) 
approximation can be used (similar to Ru- 
derman (2011a,b)). 

In this study, we are only concerned 
with the observational features of estab¬ 
lished standing slow modes, even if the 
WKB approximation is violated, i.e., the 
MHD slow mode time scale is in the same 
order of the heating or cooling time scale. 


the result could still be used to identify 
a standing slow wave based on stepwisely- 
defined quasi-equilibriums. We set up a hot 
and dense flaring coronal loop with typical param¬ 
eters that are observed by SUMER (Wang 2011). 
Our loop measures Lq = 100 Mm in length and 
a = 5 Mm in radius. The loop is filled with 
plasma with a density po = 1.4 • 10“^^kgm“^ 
(the electron density tIqq = 8.5 • 10^cm“^) and a 
temperature of T[ = 6.4 MK. The internal mag¬ 
netic field is fixed at Bi = 40 G. We choose a 
density ratio of p-Jpe = 5, a temperature ratio 
Ti/Te = 1.5, and a magnetic field strength ratio 
B^IBq = 0.91. The ratios of the plasma param¬ 
eters are set in the typical range of flaring coro¬ 
nal loops. Changes in these ratios will not sig¬ 
nificantly affect the result, since in a slow mode 
the longitudinal perturbations are strictly con¬ 
fined within the plasma cylinder and are more 
than four orders of magnitude stronger than 
the perturbations to the ambient plasma. 
In this setup, the internal and external plasma are 
typical coronal fluid with plasma beta of = 0.23 
and /3q = 0.02.6, respectively. The acoustic speeds 
are C^i = 420kms“^, Cse = 340kms“^, while the 
Alfven speed are Cai = 950kms“^ and Cag = 
2300 km s“^ (Table 1). These speeds are typical 
values observed in the solar corona (see, e.g., As- 
chwanden 2005). We also investigate the slow 
modes in hot loops at = 8.8,12,15, 20 MK 
(see Table 1), and explore how the prop¬ 
erties would deviate from the case of the 
6.4 MK loop. The total pressure pertur¬ 
bation is kept unchanged, the amplitude 
of perturbed density, temperature and ve¬ 
locity will vary according to the equillib- 
rium temperature (Table 1). However, the 
ampltiude of the perturbed quantities will 
not affect overall result, since in linearized 
MHD wave modes they are scalable. In 
the following text, we refer to the case of 
slow wave in the 6.4 MK loop by default, and 
other cases are specified as otherwise in Sec¬ 
tion 4.3. 

For n = 1 (the fundamental longitudinal mode), 
ka = 0.157 is in the long wavelength limit. We 
use A[ = 0.01 so that the velocity perturba¬ 
tion is about 57kms“^ and the density pertur¬ 
bation is about 12% of the equilibrium value. Eor 
n = 2(3), we use Ai = 0.02(0.04). The velocity 
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and density perturbation are about 41(48) kms ^ 
and 9(11)% of the equilibrium value, respectively 
(Table 1). These amplitude values are chosen 
to agree with the observed Doppler shift veloci¬ 
ties in Wang (2011). By solving the dispersion 
relation Equation 10, we obtained the periods 
Po = 520 s, 260 s, 170 s for the n = 1,2,3 modes, 
respectively. 

Starting from the equilibrium loop model, we 
construct a discrete standing slow wave model as 
specified by Equation 11-15. The simulation do¬ 
main ranges from [0,1.5a], [0,27r], [0,Po] foi* 

0-, and ^-axes, with grid cells of 150 x 180 x 300, 
respectively. Eigure 1 illustrates a snapshot of the 
p, T and Vz distribution at t = Po/8 for the n = 1 
mode. An anti-node in terms of density perturba¬ 
tion is present at z = Po/2, this is in agreement 
with SUMER observations that Doppler shift os¬ 
cillations are usually effectively detected at the 
loop apex rather than the loop footpoints (Wang 
et al. 2007). The density and temperature per¬ 
turbations are in phase and are a quarter pe¬ 
riod out of phase, with respect to both time and 
space, with the longitudinal velocity. We note 
that the density/temperature and velocity 
perturbations are in phase for propagating 
slow waves (Sakurai et al. 2002; Wang et al. 
2009b), therefore, a mix of propagating and 
standing wave would lead to the detection 
of a rather complex phase lag (Wang et al. 
2009a). In this study, we focus on purely 
standing slow MHD waves to obtain guide¬ 
lines for observations. 

3. Methods 

3.1. Spectroscopic Modelling 

We are concerned with synthesizing the ex¬ 
treme ultraviolet (EUV) emission intensity ^Ao [ergs cm ^ 
of a specific spectral line Aq for optically thin 
plasma along the LOS (Dere et al. 1997), 

ho = ^ J Gx,ine,T)nldl, (16) 

where is the abundance of the emitting element 
relative to hydrogen, Gaq [ergscm^ s“^] is the con¬ 
tribution function that contains the terms relative 
to atomic physics, dl is an infinitesimal element 
length along the LOS. 



Eig. 1.— A snapshot of p(r, z) (a), T(r, z) (b), and 
Vz{r^ z) (c) at t = Po/8. The dashed lines indicate 
the position of the rays at LOS of 30°, 45°, 60°, 
and 90°, respectively. 


To calculate the integrated emission intensity 
in Equation 16, we generate a look-up table for 
the Ee XIX AIII 8 . 1 A line using the CHIANTI 
v7.1 atomic database (Dere et al. 1997; Landi 
et al. 2013). The look-up table is sampled at 
a uniform mesh of size 200 x 200 grid points at 
logne[cm“^] G [8,11] and logT[K] G [6.2,7.7]. 
We used the CHIANTI collisional ionization equi¬ 
librium file chianti-v7.ioneq and coronal abun¬ 
dance sun-Coronal-2012-schmelz.abund (also see 
Antolin & Van Doorsselaere 2013)^. The choices 
of ionisation and abundance files do not affect our 
gre^hsjat all, since we are only concerned with the 
relative intensity perturbation and Doppler shift 
caused by the MHD waves. 

Figure 2 illustrates the forward mod¬ 
elling method and how the LOS integration 
(Equation 16) is implemented numerically. 
Eor each grid point, the emissivity = Gx^n^ 
is calculated and is spread to a Gaussian spec¬ 
trum with the width determined by the thermal 


^The source of the forward modelling code (FoMo) is avail¬ 
able at https://wiki.esat.kuleuven.be/FoMo 
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Table 1: Parameters of the loop systems and the standing slow modes 


Loops 

Ti = 6.4MK : 

li = 8.8 MK 

Ti = 12 MK 1. 

"i = 15 MK : 

li = 20 MK 

I/O [Mm] 

100 

100 

100 

100 

100 

a [Mm] 

5.0 

5.0 

5.0 

5.0 

5.0 

A[G] 

40.0 

40.0 

40.0 

40.0 

40.0 

ae[G] 

43.8 

45.3 

47.0 

48.7 

51.6 

Pi [10“^^kgm“^] 

1.4 

1.4 

1.4 

1.4 

1.4 

riei [10®cm“^] 

8.5 

8.5 

8.5 

8.5 

8.5 

PilPe 

5.0 

5.0 

5.0 

5.0 

5.0 

Ti [10®K] 

6.4 

00 

bo 

12.0 

15.0 

20.0 

Ti/Te 

1.5 

1.5 

1.5 

1.5 

2.0 

A 

0.23 

0.32 

0.44 

0.55 

0.74 

/?e 

0.026 

0.034 

0.043 

0.050 

0.044 

Csi [kms 1] 

420 

490 

570 

640 

740 

Vai [kms"i] 

950 

950 

950 

950 

950 

Cxi [kms“^] 

380 

440 

490 

530 

580 

Cse [kins“^] 

340 

400 

470 

520 

530 

Vas [kms“i] 

2300 

2400 

2500 

2600 

2700 

Cxe [kms“^] 

340 

400 

n 

460 

= 1, A = 0.01 

510 

520 

[kms“^] ^ 

57 

47 

38 

31 

28 

Pi/Pi 

0.12 

0.08 

0.06 

0.04 

0.03 

Tf/Ti 

0.08 

0.06 

n 

0.04 

= 2,A = 0.02 

0.03 

0.02 

[kms“^] 

41 





pVpi 

0.09 





T°/Ti 

0.06 

n 

= 3, A = 0.04 



[kms“^] 

48 





pVpi 

0.11 





T?/Ti 

0.07 






“The superscript 0 here and thereafter denotes the amplitude of the perturbed quantities in Equation 11-15 excluding the spatial 
and temporal terms. 
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Fig. 2.— Schematic diagram illustrates the 
forward modelling method and LOS inte¬ 
grations. 

broadening of the spectral line. The spectrum cov¬ 
ered in this study is centered at A1118.1 A with a 
range of ±0.5 A (±130km/s). This is a Fe XIX 
EUV line, in which most spectroscopic observa¬ 
tions on standing slow wave were performed (see 
Wang 2011). We used 60 wavelength values to 
sample the spectrum. It corresponds to a spectral 
resolution of AA = 16.7 mA or Av^ = 4.5km/s, 
which is sufficient to resolve the spectrum. Then 
along the LOS, the emission wavelength of the el¬ 
ementary plasma fluid is modified by the velocity 
perturbation caused by the wave and is re-binned 
into the discrete spectrum. By assuming a Gaus¬ 
sian distribution for the integrated spectrum, we 
obtain the emission intensity Ix^ and Doppler shift 
velocity vd. 

We perform the calculations for LOS angle 
0 = 30°, 45°, 60°, and 90°, respectively. The pro¬ 
jected plane-of-sky has a mesh grid of x N, = 
45 X 300 sin so that the synthesized LOS emis¬ 
sion plane has a pixel size of 0.33 X 0.33 Mm^ 
(Figure 2). 

3.2. Imaging Modelling 

To synthesize the observational features of 
SDO/AIA channels, we calculated the AIA tem¬ 


perature response function (rig, T) [ DN cm^ s ^ ] 

for bandpass a (Boerner et al. 2012): 

POO 

K^{ne^T)= / G(A,ne,T)i?,(A)dA, (17) 

Jo 

where (A) [ DN cm^ sr photon”^] is the instrument- 
wavelength response function. This is the prod¬ 
uct of the solid angle occupied by a unit sur¬ 
face (0.6" X 0.6") relative to the telescope and 
the value calculated with the aia_get_resp.pro. 
G(A, rie, T)[photoncm^ s“^ sr“^] is the contribu¬ 
tion function calculated with the isothermal.pro 
routine in CHIANTI (see Del Zanna et al. 2011). 
Then the flux Fb(x)[DNs“^] at pixel x is inte¬ 
grated along the LOS, 

F„(x) = J^Ka{ne,T)nldl. (18) 

Equation 18 is solved using the same algo¬ 
rithm as Equation 16 ( Figure 2), however, 
the intensity is obtained by summing up the 
contributing emissions in all effective wave¬ 
lengths (Equation 17) rather than spread¬ 
ing into a spectrum (see Section 3.1). A 
look-up table for each AIA bandpass is sampled 
at a uniform mesh of size 200 x 200 grid points at 
logUe [cm“^] G [8,11] and logT [K] G [4,8]. 

We synthesized the AIA 94A channel that 
would image our flaring loops (> 6.4MK). The 
resultant pixel size is kept uniform at 0.6" x 0.6", 
therefore we choose a mesh grid of Nx x N, = 

35 X 230 sin ^ for output and perform the calcu¬ 
lations for LOS angles 0 = 30°, 45°, 60° and 90°, 
respectively. Point spread functions (PSF, see 
Poduval et al. 2013; Antolin & Van Doorsselaere 
2013) would only have marginal effect at the edges 
and the cylinder boundaries. Moreover, we use 
long wavelength limits and the plasma motions 
are predominantly longitudinal, thus PSF effect is 
neglected in this study. 

4. Results 

In this section, we divide the results into two 
categories that observers mainly use to detect 
standing slow mode in a coronal loop. As we 
use linearization of the MHD wave model and 
perform geometric integration by modelling the 
atomic emissions, the results could be scaled to 
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the range of interests. In the following subsec¬ 
tions, we present a typical observation mode of 
a spectrograph or imager and study other effects 
that may affect the observations. 

4.1. Spectroscopic observations 

4 . 1 . 1 . Typical observation 

We first mimic a sit-and-stare campaign of a 
spectrograph, e.g., SUMER. The slit is ideally 
placed to fully cover the central line of our loop 
with a viewing angle of 45°. Figure 3 presents 
typical observables of a spectrograph. The inten¬ 
sity perturbation (Figure 3(a)) has a larger value 
at the loop footpoints than at the loop apex, while 
the Doppler shift vd (Figure 3(c)) shows an oppo¬ 
site spatial pattern: it has a maximum at the loop 
apex. This is in agreement with the standing slow 
wave model (see Section 2.2), in which the lon¬ 
gitudinal velocity and density are phase shifted 
in space by a quarter wavelength. Time series 
(Figure 3(b, d)) taken at a position off the loop 
apex and footpoints show a typical observation, 
which could be directly compared with Figure 2 in 
Wang (2011): the intensity perturbation oscillates 
with a quarter-period out of phase in time with 
the Doppler shift v^. We also show that intensity 
variations are in phase with line width w oscilla¬ 
tions (Figure 3(b,f)). Figure 3(e) illustrates a typ¬ 
ical spectrum that would be observed in a sit-and- 
stare mode: the spectral line is Doppler shifted 
by the standing slow wave and is also broadened 
by the perturbed velocity and temperature along 
LOS. Spectral observations at a slit position could 
provide the first signal of a MHD wave. The line 
width variation is not reported so far by any obser¬ 
vations, its temporal variation and phase relations 
with other observables could be used to identify 
the wave mode. This may be due to the low am¬ 
plitude of line width oscillation (about Ikm/s), 
thus it is beyond the detection capability of cur¬ 
rent instruments. 

When the viewing angle is normal to the loop, 
i.e., 90° (Figure 4), the intensity and line width 
broadening modulation are still significant and are 
in phase with each other. However, the Doppler 
shift vd oscillation becomes below noise level and 
is not detectable. If a slit is placed over the loop 
apex (anti-node), then the intensity and line width 
modulation will be very small. 


A measure of goodness-of-fit to a Gaus¬ 
sian spectrum is not investigated in this study 
(see e.g., Antolin & Van Doorsselaere (2013)). It 
measures the level of goodness that a combina¬ 
tion of multiple temperature spectral components 
could approximate a single Gaussian spectrum. 
For the standing slow mode, the plasma motions 
are predominately longitudinal, so plasma advec- 
tions across the cylinder boundary are negligible. 
Therefore, the spectra rarely deviate from a Gaus¬ 
sian shape. In our study, the measures at the 
order of 10“^. It means the Gaussian profiles are 
in accord with the spectra or the error variance has 
been overestimated. The latter is true in our case, 
since we do not have error bars associated with 
the spectra and unity is used as error variance. 

4 . 1 . 2 . LOS effect 

Figure 5 presents the two extreme cases at the 
loop apex and footpoints^. At the footpoints, the 
intensity and line width variations do not change 
significantly with varying LOS angles, this is be¬ 
cause the wavelength of n = 1 standing slow 
mode is significantly longer than the loop radius. 
Thus, rays of any LOS angle are less likely to tra¬ 
verse both the node and anti-node. This is con¬ 
trary to the short wavelength case for the fast 
sausage mode (Antolin & Van Doorsselaere 2013). 
In the short wavelength limit, LOS rays would 
traverse more fine structures along both radial 
and longitudinal directions, therefore the intensity 
and line width modulations are more complicated. 
However, in our case, the modulation of Doppler 
shift is significantly affected by LOS angles, since 
the longitudinal velocity dominates the perturba¬ 
tion. In a rough approximation, the amplitude 
of Doppler shift could be linearly de-projected by 
a factor of 1 — ^/90° or 1 — 20 /it , rather than 
cos^ (Figure 6). A cos^ de-projection could over¬ 
estimate the velocity perturbation by a factor of 
more than 10%. This empirical formula could be 
used in observations. The reason for the devia¬ 
tion from simple geometric projection is that each 
fluid element is projected by cos 0 , while the over¬ 
all Doppler-shifted spectrum including the contri¬ 
bution of all fluid elements along LOS does not 


^we truncated 2a/ sin 30° = 4a off the loop ends, where rays 
with LOS angle of 30° would traverse only part of the loop 
cross-section and will contain an edge effect; therefore, we 
refer to the positions of 2 ; = 4a, Lq — 4a as footpoints. 
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Fig. 3.— (a) The baseline ratio time-distance plot of the relative intensity variation {Ixq — Io)/Io along the 
central axis of the loop, (c) The time-distance plot of Vz along the central axis of the loop. The solid line 
denotes the position where we took time series of —/o)//o as shown in panel (b), (d) shows the variation 
of the Doppler shifts, (f) shows the variation of line widths and (e) shows the variation of the spectral line. 
All the information is extracted for mode n = 1 at a viewing angle 0 = 45°. 


necessarily follow the same trend. This result will 
affect the estimate of realistic longitudinal veloc¬ 
ity perturbation with a two dimensional observa¬ 
tion and the associated wave energy budget. The 
LOS effect may also lead to incorrect longitudinal 
mode identification (see, e.g., Antolin & Verwichte 
2011 ). 

The Doppler shift modulation exhibits 
similar patterns due to LOS variation at 
both the apex and footpoints (Figure 5(b, 
e)). Intriguingly, we found that positive 
modulation of the intensity is overwhelm¬ 
ingly in excess over negative modulations 
at every LOS angle (Figure 5(a) and Fig¬ 
ure 3 and 4). This effect may halve the 
periodicity at regions where the intensity 
modulation is relatively weak, e.g., at the 
apex (Figure 4(d)). This is not introduced 
by asymmetry in the geometry nor the distribu¬ 
tion of electron density Ue or plasma tempera¬ 
ture T. Figure 7 presents the contribution func¬ 
tion of All 18.1 A, the plus sign indicates the loop 
parameters. Cxoine^T) varies less than 1 % with 
the electron density Ue at high temperatures, but 
strongly peaks in temperature. In our case, pos¬ 


itive temperature modulation leads to larger in¬ 
crease in Gxq than the same amount of negative 
temperature modulation would do, therefore we 
have excess intensity enhancement when temper¬ 
ature increases. This was also found, although 
not mentioned, in the case of intensity variation 
in 193A bandpass (Antolin & Van Doorsselaere 
2013). In contrast, this effect is missing in the 
171A bandpass, as the contribution function G 171 
has almost equal gradient with respect to T in the 
temperature of interest, see Figure 6 in Antolin 
& Van Doorsselaere (2013). This effect leads to a 
halving of periodicity in intensity and line width, 
especially in small LOS angles. It would become 
more significant for a loop with a temperature such 
that dG\^ldT‘^ reach extreme values. In such con¬ 
ditions, this effect would spread to broader spatial 
regions and larger LOS angles. We also notice that 
this effect would lead the asymmetry in emission 
intensity modulation (Figure 5(a)), however, it is 
highly likely to be neglected or removed by the 
data processing technique, e.g., detrending, run¬ 
ning difference. 
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Fig. 4.— The same as Figure 3, but for n = 1, 6> = 90°. 



Fig. 6.— The Doppler shift velocity amplitude 
as a function of LOS viewing angle at the foot- 
point and apex, respectively. The dashed and solid 
lines represent the trend of cos^- and (1 — ^/90)- 
dependence, respectively. 

4 . 1 . 3 . Longitudinal overtones 

We perform sit-and-stare mode observations for 
n = 2,3 modes as well (Figures are not shown 
here). In spectroscopic observations, only a small 
segment of the spatial distribution of emission in¬ 
tensity and Doppler shift along a loop is normally 
measured. Without spatial information, it is not 
possible to judge the longitudinal mode number. 


However, with imaging observations, longitudinal 
overtones may be observed (see the next section). 

4.2. Imaging observation 

In imaging observations, the spatial distribu¬ 
tion of the emission intensity is obtained. One 
could easily follow the spatio-temporal variation of 
a loop oscillation using the time-distance method 
(see, e.g.. Yuan & Nakariakov 2012). Figure 8(a) 
shows the baseline-ratio difference plot of n = 1 
mode along the loop. The intensity perturbation 
is more significant at the loop footpoints than the 
loop apex. By changing the viewing angle, the in¬ 
tensity modulation changes slightly, but the over¬ 
all spatio-temporal pattern remains unchanged. 
This is in agreement with the spectroscopic model, 
see discussions in Section 4.1. The time-distance 
plot along a loop could reveal the nodal structures 
of the longitudinal overtones (Figure 8(b,c). By 
comparing the amplitude variation along the loop, 
one could measure the longitudinal mode number 
n, and hence the wavelength of the slow mode 
2I/o/n. 

4.3. Observing standing slow mode in hot 
loops with different temperatures 

So far, we only studied a 6.4 MK loop 
with fixed density and temperature ratio. 
As long as the background emission is in- 
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Fig. 5.— One period variation of the normalized emission intensity I/Iq (a, d), Doppler shift velocity vd 
(b, e), and line-width re (c, f) observed at the loop footpoint and apex for the n = 1 mode. 





significant compared to loop emission, the 
current result remains valid. However, if a 
coronal loop is heated to the hotter wing 
(Ti > 8.8 MK) of Gxo (Figure 7), the back¬ 
ground emission could become stronger 
than the loop itself. In such cases, the 
observed Doppler shift and intensity mod¬ 
ulation would be rather different. The free 
parameters are the loop temperature, den¬ 
sity, the temperature and density ratio. It 
is not practical to iterate all possible combi¬ 
nations, so we only present possible scenar¬ 
ios as how to understand a slow standing 
mode properly. We vary the loop temper¬ 
ature, while keeping other free parameters 
unchanged. Several loop models are set 
up at the various temperatures where the 
SUMER A1118.1A responses are signficant 
(see Table 1). 

Figure 9 illustrates snapshots of the rel¬ 
ative intensity, Doppler velocity and line 
width of standing slow modes observed in 
loops at Ti =6.4, 8.8, 12, 15, and 20 MK, 
respectively. The snapshots are taken at 
t = Pq/S at a viewing angle 0 = 45° 
(Other viewing angles give similar results, 
therefore they are not shown). At the 
cooler wing of AIII 8 .IA, positive tempera¬ 
ture modulation would enhance the inten¬ 
sity, while at the hotter wing, the oppo¬ 
site would occur. The intensity modula¬ 
tion becomes relatively small, as the am¬ 
bient plasma emission intensity become sig¬ 


nificant or stronger than the loop itself (e.g, 
Ti = 20 MK). 

Figure 10(b) shows the normalized rel¬ 
ative intensity modulation {51^/Iq)/{ 5p/po)^ 
which is usually assumed to be 2 (Wang 
et al. 2009b), at the loop footpoint against 
loop temperature. It illustrates how the 
sign of dGxo{T)/dT would modulate the in¬ 
tensity variation and causes the asymme¬ 
try effect (see discussions in Section 4.1.2). 
We found that {61^/Iq)/{ 6p/po) could be con¬ 
siderably larger than 2, and that it is only 
close to 2 at regions where \dGxo{T)/dT\ is 
relatively small. At extremely hot loops 
(20MK), the normalized intensity modula¬ 
tion could even approach zero. 

On the other hand, the Doppler shift 
velocity of slow standing mode remains de¬ 
tectable for most loop temperatures. How¬ 
ever, the normalized amplitude COS 9) 

at the loop apex (it is symmetric for both 
positive and negative motions) deviates 
more and more from unity for loops at 
higher temperatures (Figure 10(a)). This 
does not mean that the wave energy is un¬ 
detected. It is still buried in the line width 
(Figure 9, right column) or the skewness of 
the spectra, which are not measured in this 
study. 
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Fig. 7.— The contribution function Gx^iriQ^T) for the Fe XIX AIII 8 . 1 A line. The maximum formation 
temperature is logT = 6.95 (T = 8.9 MK). The plus symbol marks the region with the used loop parameters 

[riehTi]. 


5. Conclusion 

In this paper, we performed forward modelling 
to predict the observational features of standing 
slow modes in hot coronal loops. By considering 
the geometrical and instrumental effects, we mea¬ 
sure the plasma emission intensity, Doppler shift, 
line width and spectrum modification caused by 
the standing slow mode. 

We found that the Doppler shift velocity is sig¬ 
nificantly affected by LOS effect, it could become 
undetectable at an viewing angle of ^ = 90°. We 
found that a linear de-projection by the LOS an¬ 
gle is more accurate than the cos6> de-projection 
as traditionally used in Wang (2011). The emis¬ 
sion intensity perturbation is normally a quarter- 
period out of phase with the Doppler shift ve¬ 
locity variation, both in time and space. This 
effect has been used to identify the standing slow 
wave mode. Positive temperature variation in¬ 
troduces more emission intensity enhancement 
than the same amount of negative temperature 
variation would reduce it, because dCx^/dT\rj.+ 
is significantly larger than 9 Gao/ 5T|^-, i.e., 
d‘^Gxo/dT‘^\TQ reaches positive extremes. This 
effect could lead to the halving of the peri¬ 
odicity in intensity and line width at the loop 
apex. Half periodicity could be also reached if 
the OGxq/ dT\rj.- > dG Xq / dT\rj.+ , which could be 


only found at the hotter wing of Gaq, meaning 
that the loop has to be heated to a few tens of 
MK in mega flares. This second-order effect in 
the contribution function could also lead to asym¬ 
metry in the emission intensity modulation. 

Spectroscopic observations with a sit-and-stare 
mode alone are not able to resolve longitudinal 
overtones along the loop, due to the lack of spatial 
information. With EUV imagers, the longitudinal 
overtones could be resolved by investigating the 
spatial distribution of the emission intensity mod¬ 
ulation. A good way of studying standing slow 
modes is to use joint observations of spectrographs 
and imagers. 

The Doppler shift oscillation of a standing slow 
wave strongly depends on the viewing angle, for 
those observations with LOS angle close to 90°, 
we may not detect Doppler shifts. For loops on 
the solar disk, a sit-and-stare campaign of spec- 
trographic observation is favoured to be placed 
slightly off the loop apex, because there is a higher 
probability that the viewing angle would be close 
to 90°. For spectrographic observations off-limb, 
most loops apexes are well exposed for observation 
at good viewing angles. 

The emission contribution function plays 
a signficant role in determining the obser¬ 
vational features of a standing slow wave. It 
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Fig. 8.— Baseline ratio time-distance plot of AIA 94 A emission intensity along the central axis of the loop 
for modes n = 1 (a), 2 (b), 3 (c), respectively. 
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(b) (l-lo)/lo, n=2 


may cause asymmetric intensity modulation 
for positive and negative temperature per¬ 
turbations. This effect could be the oppo¬ 
site for loops at the cooler and hotter wings 
of the AIII 8 . 1 A line. The normalized rel¬ 
ative intensity modulation {61^/Iq)/{ 6p/po) 
would be 2 only for the loops close to the 
peak response temperature of the A1118.1 A 
line. The Doppler shift velocity could be 
significantly smaller than the plasma veloc¬ 
ity, if the background plasma emission be¬ 
comes more significant or stronger than the 
loop emission. For loops at the hotter wing 
of a spectral line, the intensity modulation 
could be small, and one will only observe 
the Doppler shift. It may lead to a false in¬ 
terpretation of the result, e.g., in the case 
of a propagating fast wave Tomczyk et al. 
(2007); Van Doorsselaere et al. (2008). 

Imaging observations are subject to more vari¬ 
ability, e.g., loop curvatures, plasma stratification, 
ambiguity of the loop footpoints. The amplitude 
variation along a loop is sufficient to identify the 
longitudinal overtones (Figure 8). Normally the 
n = 1 mode could be well identified even if the 
footpoints are hard to find. However, the n > 1 
overtones have shorter wavelength and the nodes 


close to the footpoints are less reliably measured, 
therefore, it poses some challenges to identify a 
n > 1 overtone. 

From the observables, we could attempt to per¬ 
form MHD seismology (Wang et al. 2007). For the 
input parameters, we allocate 10% uncertainties 
to period Pq, loop length Lq. While we keep the 
perturbations of density and temperatures as the 
uncertainties. So we have loop length Lq = 100 ± 
10 Mm, density po = (1.4±0.2) • 10“^^ kgm“^ and 
temperature Tq = (6.3±0.5) MK, then we estimate 
the magnetic field strength Bq = (40.6 ±6.2) G in¬ 
side the loop. If we compare the result to the real 
input 40 G, the main uncertainty is in the estima¬ 
tion of the parameters, while the assumption that 
2 L 0 IP 0 Ct only result in about 1.5% of uncer¬ 
tainty. In the long wavelength limit, this assump¬ 
tion remains valid and causes small uncertainties 
in MHD seismology. However, one has to be cau¬ 
tious in using this relation, if the wavelength is 
much shorter or the plasma p is not small, 
the associated relative error could be esti¬ 
mated by (1 ±7/3i/2)^-^ — 1. 

The research was supported by an Odysseus 
grant of the FWO Vlaanderen, the lAP P7/08 
GHARM (Belspo), the Topping-Up grant Gor- 
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respectively. 


Seis and the GOA-2015-014 (KU Leuven), and the 
Open Research Program KLSA201312 of Key Lab¬ 
oratory of Solar Activity of National Astronomi¬ 
cal Observatories of China (D.Y.). CHIANTI is a 
collaborative project involving Ceorge Mason Uni¬ 
versity, the University of Michigan (USA) and the 
University of Cambridge (UK). 

REFERENCES 

Andries, J., Arregui, L, & Goossens, M. 2005, 
ApJ, 624, L57 

Anfinogentov, S., Nakariakov, V. M., Math- 
ioudakis, M., Van Doorsselaere, T., & Kowal¬ 


ski, A. F. 2013, ApJ, 773, 156 

Antolin, R, & Rouppe van der Voort, L. 2012, 
ApJ, 745, 152 

Antolin, P., & Van Doorsselaere, T. 2013, A&A, 
555, A74 

Antolin, P., & Verwichte, E. 2011, ApJ, 736, 121 

Antolin, P., Vissers, G., Pereira, T. M. D., Rouppe 
van der Voort, L., & Scullion, E. 2015, ArXiv 
e-prints, arXiv: 1504.04418 

Antolin, P., Yokoyama, T., & Van Doorsselaere, 
T. 2014, ApJ, 787, L22 


15 














































0 = 45 ° 




T, [MK] 


Fig. 10.— For loops at different temperatures, (a) 
presents the normalized doppler velocity at a view¬ 
ing angle of ^ = 45° at loop apexes; (b) shows the 
normalized relative intensity modulation at loop 
footpoint, where 61°/Iq = {It=o-It=Po/ 4 ,)/It=Po/ 4 , 
and —5I°/Io = {It=Po /2 — It=Pol4:)IIt=Pol4:- The 
shade areas illustrate the scaled Gxq{T) (a) and 
dGxQ{T)/dT (b) at Ue = 8.5 • 10^cm“^, respec¬ 
tively. The diamond and asterisk denote, respec¬ 
tively, the positive and negative modulations for 
both Doppler shift and relative intensity. 


Arregui, I. 2015, ArXiv e-prints, arXiv: 1501.06708 

Arregui, L, Andries, J., Van Doorsselaere, T., 
Goossens, M., & Poedts, S. 2007, A&A, 463, 
333 

Asai, A., Shimojo, M., Isobe, H., et al. 2001, ApJ, 
562, L103 

Aschwanden, M. J. 2005, Physics of the Solar 
Corona. An Introduction with Problems and 
Solutions (2nd edition) 


Boerner, P., Edwards, C., Lemen, J., et al. 2012, 
Sol. Phys., 275, 41 

Chen, Y., Feng, S. W., Li, B., et al. 2011, ApJ, 
728, 147 

Chen, Y., Song, H. Q., Li, B., et al. 2010, ApJ, 
714, 644 

Cooper, F. C., Nakariakov, V. M., & Tsiklauri, D. 
2003a, A&A, 397, 765 

Cooper, F. C., Nakariakov, V. M., & Williams, 
D. R. 2003b, A&A, 409, 325 

De Moortel, L, Hood, A. W., Ireland, J., & Walsh, 
R. W. 2002a, Sol. Phys., 209, 89 

De Moortel, L, Ireland, J., & Walsh, R. W. 2000, 
A&A, 355, L23 

De Moortel, L, Ireland, J., Walsh, R. W., & Hood, 
A. W. 2002b, Sol. Phys., 209, 61 

De Moortel, L, & Nakariakov, V. M. 2012, Royal 
Society of London Philosophical Transactions 
Series A, 370, 3193 

Del Zanna, G., O’Dwyer, B., & Mason, H. E. 2011, 
A&A, 535, A46 

Dere, K. P., Landi, E., Mason, H. E., Monsignori 
Eossi, B. C., & Young, P. R. 1997, A&AS, 125, 
149 

Edwin, P. M., & Roberts, B. 1983, Sol. Phys., 88, 
179 

Eleishman, G. D., & Kuznetsov, A. A. 2010, ApJ, 
721, 1127 

Goossens, M., Erdelyi, R., & Ruderman, M. S. 
2011, Space Sci. Rev., 158, 289 

Gruszecki, M., Nakariakov, V. M., & Van Doors¬ 
selaere, T. 2012, A&A, 543, A12 

Hood, A. W., Browning, P. K., & van der Linden, 
R. A. M. 2009, A&A, 506, 913 


Aschwanden, M. J., Eletcher, L., Schrijver, C. J., 
& Alexander, D. 1999, ApJ, 520, 880 

Aschwanden, M. J., Nightingale, R. W., Andries, 
J., Goossens, M., & Van Doorsselaere, T. 2003, 
ApJ, 598, 1375 


Jess, D. B., Reznikova, V. E., Van Doorsselaere, 
T., Keys, P. H., & Mackay, D. H. 2013, ApJ, 
779, 168 

Kim, S., Nakariakov, V. M., & Shibasaki, K. 2012, 
ApJ, 756, L36 


16 






Klimchuk, J. A. 2006, Sol. Phys., 234, 41 

Krishna Prasad, S., Banerjee, D., & Van Doorsse- 
laere, T. 2014, ApJ, 789, 118 

Kumar, R, Innes, D. E., & Inhester, B. 2013, ApJ, 
779, L7 

Kumar, P., Nakariakov, V. M., & Cho, K.-S. 2015, 
ArXiv e-prints, arXiv:1502.07117 

Kupriyanova, E. G., Melnikov, V. E., & Shibasaki, 
K. 2013, Sol. Phys., 284, 559 

Kuznetsov, A. A., Van Doorsselaere, T., 
& Reznikova, V. E. 2015, Sol. Phys., 
arXiv:1502.06716 

Landi, E., Young, P. R., Dere, K. P., Del Zanna, 
G., & Mason, H. E. 2013, ApJ, 763, 86 

Mariska, J. T. 2006, ApJ, 639, 484 

Melnikov, V. E., Reznikova, V. E., Shibasaki, K., 
& Nakariakov, V. M. 2005, A&A, 439, 727 

Moreels, M. G., & Van Doorsselaere, T. 2013, 
A&A, 551, A137 

Mossessian, G., & Eleishman, G. D. 2012, ApJ, 
748, 140 

Nakariakov, V. M., & Melnikov, V. E. 2009, 
Space Sci. Rev., 149, 119 

Nakariakov, V. M., & Ofman, L. 2001, A&A, 372, 
L53 

Nakariakov, V. M., Ofman, L., Deluca, E. E., 
Roberts, B., & Davila, J. M. 1999, Science, 285, 
862 

Nakariakov, V. M., & Verwichte, E. 2005, Living 
Reviews in Solar Physics, 2, 3 

Ning, Z. 2014, Sol. Phys., 289, 1239 

Ofman, L. 2010, Living Reviews in Solar Physics, 
7,4 

Ofman, L., & Wang, T. 2002, ApJ, 580, L85 

Parnell, C. E., & De Moortel, 1. 2012, Royal Soci¬ 
ety of London Philosophical Transactions Series 
A, 370, 3217 

Peter, H., Bingert, S., Klimchuk, J. A., et al. 2013, 
A&A, 556, A104 


Poduval, B., DeEorest, G. E., Schmelz, J. T., & 
Pathak, S. 2013, ApJ, 765, 144 

Priest, E. R., Heyvaerts, J. E., & Title, A. M. 
2002, ApJ, 576, 533 

Reale, E. 2014, Living Reviews in Solar Physics, 
11,4 

Reznikova, V. E., Antolin, P., & Van Doorsselaere, 
T. 2014, ApJ, 785, 86 

Reznikova, V. E., Van Doorsselaere, T., & 

Kuznetsov, A. A. 2015, A&A, 575, A47 

Ruderman, M. S. 2011a, A&A, 534, A78 

—. 2011b, Sol. Phys., 271, 41 

Ruderman, M. S., & Erdelyi, R. 2009, 

Space Sci. Rev., 149, 199 

Sakurai, T., Goossens, M., & Hollweg, J. V. 1991, 
Sol. Phys., 133, 227 

Sakurai, T., Ichimoto, K., Raju, K. R, & Singh, 
J. 2002, Sol. Phys., 209, 265 

Scullion, E., Rouppe van der Voort, L., Wede- 
meyer, S., & Antolin, P. 2014, ApJ, 797, 36 

Selwa, M., Murawski, K., & Solanki, S. K. 2005, 
A&A, 436, 701 

Taroyan, Y., & Bradshaw, S. 2008, A&A, 481, 247 

Taroyan, Y., & Erdelyi, R. 2009, Space Sci. Rev., 
149, 229 

Taroyan, Y., Erdelyi, R., Doyle, J. G., & Brad¬ 
shaw, S. J. 2005, A&A, 438, 713 

Taroyan, Y., Erdelyi, R., Wang, T. J., & Brad¬ 
shaw, S. J. 2007, ApJ, 659, L173 

Tomczyk, S., McIntosh, S. W., Keil, S. L., et al. 
2007, Science, 317, 1192 

van der Holst, B., Sokolov, 1. V., Meng, X., et al. 
2014, ApJ, 782, 81 

Van Doorsselaere, T., De Groof, A., Zender, J., 
Berghmans, D., & Goossens, M. 2011a, ApJ, 
740, 90 

Van Doorsselaere, T., Gijsen, S. E., Andries, J., & 
Verth, G. 2014, ApJ, 795, 18 


17 



Van Doorsselaere, T., Nakariakov, V. M., & Ver- 
wichte, E. 2008, ApJ, 676, L73 

Van Doorsselaere, T., Wardle, N., Del Zanna, G., 
et al. 2011b, ApJ, 727, L32 

Verth, G., Erdelyi, R., & Jess, D. B. 2008, ApJ, 
687, L45 

Wang, T. 2011, Space Sci. Rev., 158, 397 

Wang, T., Innes, D. E., & Qiu, J. 2007, ApJ, 656, 
598 

Wang, T., Solanki, S. K., Gurdt, W., Innes, D. E., 
& Dammasch, 1. E. 2002, ApJ, 574, LlOl 

Wang, T. J., Ofman, L., & Davila, J. M. 2009a, 
ApJ, 696, 1448 

Wang, T. J., Ofman, L., Davila, J. M., & Mariska, 
J. T. 2009b, A&A, 503, L25 

Wang, T. J., Solanki, S. K., Gurdt, W., et al. 
2003a, A&A, 406, 1105 

Wang, T. J., Solanki, S. K., Innes, D. E., & Gurdt, 
W. 2005, A&A, 435, 753 

Wang, T. J., Solanki, S. K., Innes, D. E., Gurdt, 
W., & Marsch, E. 2003b, A&A, 402, L17 

Williams, D. R., Mathioudakis, M., Gallagher, 
R T., et al. 2002, MNRAS, 336, 747 

Williams, D. R., Phillips, K. J. H., Rudawy, P., 
et al. 2001, MNRAS, 326, 428 

Winebarger, A. R., Warren, H., van Ballegooijen, 

A. , DeLuca, E. E., & Golub, L. 2002, ApJ, 567, 
L89 

Yuan, D., & Nakariakov, V. M. 2012, A&A, 543, 
A9 

Yuan, D., Nakariakov, V. M., Huang, Z., et al. 
2014a, ApJ, 792, 41 

Yuan, D., Pascoe, D. J., Nakariakov, V. M., Li, 

B. , & Keppens, R. 2015, ApJ, 799, 221 

Yuan, D., Sych, R., Reznikova, V. E., & Nakari¬ 
akov, V. M. 2014b, A&A, 561, A19 


This 2-column preprint was prepared with the AAS 
macros v5.2. 


18 



